Temperature dependence of the nodal Fermi velocity in layered cuprates 
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We explain recently observed linear temperature dependence of the nodal Fermi velocity vp (T) 
in near-optimally doped cuprates. We argue that it originates from electron-electron interaction, 
and is a fundamental property of an arbitrary 2D Fermi liquid. We consider a spin-fermion model 
with the same parameters as in earlier studies, and show that the T term is about 30% at SOOA', 
in agreement with the data. We show that the sub-leading term in vp{T) is a regular (and small) 
T^ correction. We also show that at a 2kF quantum-critical point, temperature corrections to the 
dispersion are singular. 
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The origin of strong deviations from the Fermi liquid 
behavior in the normal state of the hole-doped cuprates, 
the mechanism of d— wave superconductivity, and the na- 
ture of the pseudogap phase remain the subjects of active 
debate in the condensed-matter community. Deep inside 
the pseudogap phase the cuprates are Mott-Hubbard in- 
sulators. Outside the pseudogap phase Angle-Resolved 
Photoemission Spectroscopy (ARPES) and other mea- 
surements show a large, Luttinger Fermi surface, and w^ 
behavior of the fermionic damping at the lowest ener- 
gies [l|, y , consistent with the idea that in this range the 
system is a Fermi liquid with strong correlations. 

We take a point of view that the crossover from a metal 
to a Mott insulator occurs inside the pseudogap phase, 
while to the right of the T* line, the number of carriers is 
1—x, where x is doping. In the 1—x regime, the fermionic 
self-energy can be described in conventional terms, as 
originating from the interaction with some bosonic de- 
grees of freedom. A boson can be a phonon, or it can be 
a collective electronic excitation in spin or charge chan- 
nels. The same interaction is also thought to be primarily 
responsible for the pairing instability, which eventually 
leads to a superconductivity. 

The nature of the pairing boson in the cuprates is 
the subject of outgoing debate, and in recent years 
several proposals to distinguish experimentally between 
phononic and electronic mechnisms have been dis- 
cussed [3|, y, la, la, lZI, y . One of the proposals is to look 
at the temperature dependence of the the Fermi velocity 
vf{T) = vf{T = 0)(1 -I- 5{T)), taken along the diago- 
nal of the Brillioine Zone (BZ) where the d^2_j,2— wave 
superconducting gap has nodes (the "nodal" velocity) 
It has been measured recently by Plumb et al. [9] in 
optimally-doped Bi2Sr2CaCu2 08+5 be means of laser- 
based angle-resolved photoemission spectroscopy. They 
found that 5{T) is approximately linear in T up to at 
least ZQQK . The linear behavior holds down to Tc, and 
the slope is quite large: 5{T) is about 0.35 between Tc 



and T = 250K. The hnear T dependence of S{T) is a 
challenge to theorists, as on general grounds one would 
expect an analytic, T^ dependence at the lowest T. The 
magnitude of 6{T) is another challenge, as the coupling 
to the boson is set by fits to other experimental data, 
including the value of Tc- 

The linear in T dependences of various observables in 
the normal (non-pseudogap) state have been reported 
before and phenomenologically described in terms of 
marginal Fermi liquid (MFL) behavior [l^, [ll[. The 
linear in T behavior of the velocity is different by two 
reasons. First, the self-energy which yields S{T) ex T 
scales as ReT, ex ujT, while MFL behavior originates from 
i?eE ex w log max (w,T) (and IrnE ex max {uj,T)). Sec- 
ond, the linear behavior of S{T) has been measured down 
to energies (ttT —f to) where ImT, already has a Fermi 
liquid, Lo^ form [9|. From this perspective, the linear in 
T dependence of the velocity appears to be a truly low- 
temperature asymptotic behavior, as opposed to MFL 
behavior, which likely holds only above the upper bound- 
ary of the Fermi liquid [Tl, [l^] . 

The quasi-linear T dependence of the Fermi velocity 
can be obtained due to electron-phonon iteraction, but 
only as approximate behavior in a limited T range above 
the Debye frequency lud UM^ where the temperature- 
dependent velocity slowly approaches its bare value af- 
ter passing through a minumum at T ~ tj/j. On the 
other hand, electron-electron interaction in 2D gives rise 
to a linear in T dependence of S{T) down to the lowest 
temperatures, at which the system is still in the nor- 
mal state la, ll6|. The linear in T correction to the 



velocity comes from cuT term in the real part of self- 
energy, which in turn originates from the non- analytic 
ImT,{uj) ~ uj^loguj, coming from backscattering lj,|l5|. 



Although the original analysis in [15], [16| was a weak- 
coupling perturbation theory, one can show that S{T) is 
linear in T in an arbitrary Fermi liquid. The reasoning 
parallels the one in Ref. [ITI for the specific heat coeffi- 



cient, which is also hnear in T in a 2D Fermi hquid. 

The magnitude of 5T is a different issue. To second or- 
der in the interaction U{q), the velocity renormalization 
is given by [l5( 



6{T) = AT 



V(0) - lu{2kp) 




(1) 
where A = pp\og2/{ATT^Vp). The two terms in 
Eq.([T]) account for the contributions from charge and 
spin channels. Approximating the measured FS of 
Bi2Sr2CaCu208+A (Bi2212) by a circle with pF ~ V2 x 
O.Gtt — 2.7 (Ref. [iSJ, we set interatomic distance to one), 
U{q) by [/ ~ 2eV, and using experimental vp ~ leF 
(Ref. [190, we find S{T) - 1.610-^r, where T is measured 
in Kelvins. This yields 5{T) - 5 x 10"^ for T - 300A^ 
two orders of magnitude smaller than the experimental 
value. 

The second-order estimate, however, is only valid in 
the weak coupling regime, when the system is far from 
a Pomeranchuk-type instability. There is a consensus 
among researchers that near optimal hole doping, the 
coupling is actually strong. Strong interaction brings 
about two effects. On one hand, larger coupling leads 
to a larger S{T). On the other hand, the effective, 
screened interaction between quasiparticles deviates from 
the second-order U^H form, where 11 is the dynamic po- 
larization bubble, and one has to consider the full effect 
of the screening. Since H is temperature-dependent, this 
generaly affects the functional form of S{T). The pur- 
pose of this paper is to analyze the interplay between 
these two effects. 

As a model for the dynamic screening, we consider 
the interaction between fermions and their collective 
bosonic excitations in the spin channel. Several groups 
argued 12, l20|j l21| that this interaction is responsible for 
the normal state self-energy and the pairing. We show 
that the velocity renormalization is not affected by the 
full dynamic screening of the effective interaction. The 
extra terms in 6(T) due to the difference between the 
full interaction and U^H, scale as T^ and are small for 
all experimentally studied T. For the linear in T term, 
we find, using the same parameters as in earlier studies, 
S{T) ss 0.37 at 250K, in a very good agreement with the 
data. 

The singular interaction in the spin channel is a 2kF 
process, which for a hole-like cuprate Fermi surface is 
an umklapp scattering between nodal fermions at kp = 
(/cfj^f) and Q — hp- The normal state fermionic self- 
energy due to nodal-nodal interaction has been previ- 
ously considered in the context of quantum-critical (QC) 
phenomena [22, l23| . Here we consider the system behav- 
ior away from the QC point, when the bosonic propa- 
gator has a finite mass. In the notation G~^ (k, iu>m) = 
itOrn — fk+kp + ^O^i ^m), the sclf-eucrgy due to coupling 



to spin fiuctuations is given by 
n 



E(k,iw™) = -^T^ /d\x(q-2kF,ifi™) 



1 



«(w,„ + f^m) — S-kF+k+q+Q 



(2) 



where kp is the Fermi momentum along zone diago- 
nal, counted from (tt, tt) point [(O.Svr, O.Stt) for optimally 
doped Bi2212]. Near the Fermi surface 

£k+kF = vpk^ + (3'^ky, 
e-kp+k+q+Q = -vp{kx + q^) + 0^{ky + qyY (3) 

where x is set along the zone diagonal, towards (— tt, — tt), 
and (3 parametrizes the curvature of the Fermi line, 0^ — 
\/{2m) for a circular FS. 

Further, x(qj^m)j normalized to x(0,0) = 1, is the 
dimensionless dynamical spin susceptibility, and g — 
{U{2kp)/2) * K is the effective, enhanced coupling in the 
spin channel. The factor K, which reduces to 1 at weak 
coupling, is the ratio of the actual and bare static spin 
susceptibilities at momentum 2kp. We assume, based on 
neutron scattering data, that the momentum dependence 
of the static susceptibility is weak near q = 0, and ap- 
proximate the full dynamic spin susceptibility at small q 
as x(q- 2'kp,i^rn) = 1/(1 -|- 5n(q - 2kp,ii}m)), where 
n(q— 2ki?, if^m) is the polarization bubble with momenta 
near 2ki? (the spin factor of 2 is included into 11) . 

The 2fc p p article-hole bubble has been calculated 
before [22I |23J . It contains a regular part, which plays 
no role in our analysis and which we neglect, and a non- 
analytic part (a dynamic Kohn anomaly), which at finite 
T is given by 



nAfA(q- 2kp,inr. 



1 



AttvpP 
At T = 0, this reduces to 



du^ ^—5 — (4) 



4Tcosh- 



2 M--Eq 

4T 



U'^-^'{ci-2kp,in,r^) 



2^^^!^'^ 



^l+El (5) 



where E^ = -I'Ffe + /3^g.y/2 [24]. The static Hlff{q - 
2ki?,0) is nonzero at the smallest q only for q^, < (a 
static Kohn anomaly). The non-analycity that gives rise 
to our effect originates from the form of n|)^° at Eq < 0, 
and fl'^ <$: E^. In this fimit Il%=f ex |^m|/v^|Sq|, which 
for qrc ~ qy is of order |ri„i|/|(7.y|. The existence of non- 
analytic \qy\ in the denominator implies that ^Jj^ gives 
rise to dynamic long-range interaction between fermions. 
For the velocity renormalization, we need the real part 
of the self-energy ReS(k, cu) on the mass shell, to first or- 
der in fermionic frequency uj: Ke'E^Sk = uj,uj) = wA(T). 
The temperature variation of the velocity is related to 



A(r) as 



the Fermi surface, and is not sensitive to a cutoff. 



5{T) 



vf{T)-vf{T = 0) A(T)-A(0) 



vf{T = 0) 



1 + A(0) 



(6) 



where vpiT = 0) = vp/{l + X{0)). The zero-temperature 
A(0) is cutoff-dependent [23|. The temperature depen- 
dence of A on the other hand comes from processes near 



We first obtain fmS](k, lj) using spectral representa- 
tion, and then obtain ReS by Kramcrs-Kronig transfor- 
mation. Substituting Eqs. ©-(H]) into Eq. ([2]), using 
spectral representation, evaluating the frequency sums, 
and re-scaling, we get 



6 T^ f^°° r 

ImS(k,w) = -^ / dt[nB{t) +nF(t + Ld)] 



l/Vt 



+ OC 



dz 



sh^z 



dylmZ 



1/2 



RcZ ) 



^^"\^z\ 



Az + t + LO + k- 2/2/2 +J{Az + t + uj + k- y^/2f - [t + i5Y , 



(7) 



(8) 



where t = Q/T is the running dimensionless frequency 
variable, y oc qyjyT is the dimensional momentum vari- 
able along the Fermi surface, np.Bit) = (e* ± 1)~^ are 
Bose and Fermi functions, respectively. L is the upper 
limit of the momentum integral along the Fermi surface, 
Ti = {A-KVFP/gf, and w = w/T, k = ypk/T. The 
second-order perturbative result is obtained by neglect- 
ing Z in denominator of ([7])- 

If the integral over y was convergent, the result of 
y— integration would be 0{t). ImS would then be de- 
termined by t = 0(1) and have a Fermi liquid form, 
w^ + (ttT)'^. ReE would then only contain a regular 
ujT'^ term, and the renormalization of the Fermi velocity 
would be T^. On a more careful look, however, we find 
that at large y, Z scales as —2it/y, and the y integral 
is logarithmic. As t/y oc \Omega\/\qy\, this logarithmic 
singularity indeed originates from long-range dynamical 
interaction given by Una- 

We assume and then verify that the logarithmic ac- 
curacy is sufficient for the T-dependence of the veloc- 
ity renormalization, and that the logarithm is cut at the 
lower end by \t\. To this accuracy, ImE(k, cj) is indepen- 
dent on fc, even in frequency, and the frequency depen- 
dence conies only via the Fermi function in ([7]) . Evaluat- 
ing the integral over y in ^ with logarithmical accuracy 
and subtracting regular lu'^ + i'^T)'^ terms, we obtain for 
the non-analytic part of ImE 



at small lu/T 



dtx 



37^2 rL^/T 

lTnT,{u;)NA = -;^ / 

t log t^ [npit + to) + npit - Co) - 2nF (t)] (9) 

Substituting ^ into Kramers-Kronig (KK) relation 
ReS((I') = (2(D/7r) /(j°° Iml](s)/(s2 — di^) and integrating 
over s, we obtain after straightforward calculations that 



ReE(w) 



6ujT 



L^ IT 



tanh tdt 



6£2 

Ti 



6 log 2 

Ti 



ujT 



(10) 
The first term is a cutoff-dependent zero-temperature 
contribution A(0)w. The second term, on the other hand, 
is a universal, cutoff-independent loT term which gives 
rise to a linear in T correction to the velocity: 



b{T) = 



6 log 2 T 



1 + A(0) Ti 



(11) 



We emphasize that this universal term conies from the 
upper limit of the integration over frequency variable t, 
justifying our assumption that the logarithmic integral 
over y is cut by t rather than by external parameters (D 
and vpk. 

To estimate the slope, we recall that g = ([/(2fc_F)/2) * 
K, where K is the ratio of the actual and bare static 
spin susceptibilities at momentum 2kp. The bare sus- 
ceptibility xo = Pf/t^vf ~ 0.9states/eV. The mea- 
sured X lias a flat top between (tTjTt) and 2kF, at about 
ISstates / eV [25[. Using the ratio as the estimate for K, 
and U{2kF) ~ 2eV, pp '^ \/2 * 0.67r/a = 2.7, vp ~ leV 
[Ep - 1.35el/), we obtain Ti - 2.5 * WK (Ref. |2^. 
Using next A(0) = 0.7 extracted from ARPES fits [27|, 
we find 5T « 0.15 x lO'^T. For T = 250K, this gives 
S{T) « 0.37, in a good agreement with the data. 

We next consider the subleading terms. By power 
counting, the subleading terms, obtained by expanding 
in Z from the denominator in Eq.©, should scale as 
\JT jTx . As T\ ^ 2500i4r and the measurements are 
performed up to T ~ SOOiiT, such corrections would be 
substantial. We argue, however, that the corrections to 
Eg. pT)) are in fact regular T^ terms. To see this, we ob- 
serve that both ReZ and Im^ in Eq. (O vanish at large y. 
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FIG. 1: (color online) (a) The real part of the self-energy 
for various temperatures in the normal state (from the nu- 
merical solution of Eqs.((7])-(l8](l7]) and Kramers-Kronig trans- 
form) . The inset shows the imaginary part of the self-energy, 
(b) The temperature dependence of the coupling constant 
A(T) = dRe'£{uJ -^ 0,T)/cLj. The straight line is the ana- 
lytical result at low T, Eq. [TT|) . The inset shows the A(r) for 
three different L. 



such that expanding the denominator in ^ in powers of 
Z and integrating over y, wc loose the logt^ term, which 
was the source of non-analyticity. Integrating further 
over t we find that the expansion of l'niE{uj) is regular 
and holds in even powers of cj^ and T^, in which case the 
renormalization of the Fermi velocity holds in powers of 
T^. Alternatively speaking, the T term in vf{T) comes 
from ImZ term in Eq.([7]), while Z-dependent terms in 
the denominator in ([7]) only give rise to regular, T^ cor- 
rections to the Fermi velocity. 

To estimate these regular corrections, we computed 
the velocity renormalization numerically. We used Eq. 
[7] as a point of departure, evaluated Iml](fc,cj) by ex- 
plicit 3D integration, and subtracted antisymmetric con- 
tributions to ImS which do not affect linear in lo term 
in ReE (such contributions are generally present because 
the dispersion is not particle- hole symmetric). In Fig. [1^ 
we present the results for the real part of the self-energy 
for various temperatures, measured in units Ti. We see 
that the slope of the real part of the self-energy decreases 
with increasing temperature and is negative. In Fig. [T)d 
we plot the T-dependent part of A(r) vs Eq. (HI]). We 
see that the agreement is perfect up to approximately 
O.lTi. At larger T, the actual A(T) flatters. In the insert 
to this figure, we show that the T dependent part of A 
(but not A(0)) is unsensitive to the upper cutoff of the 
momentum integration, in agreement with Eq. (Ilip . 

The linear in T dependence of the Fermi velocity 
cannot be carried over to the 2kF QC, non-Fermi liq- 
uid regime, because the pre-factor for the T term con- 
tains the divergent static spin susceptibility at 2kp (via 
g oc K). This quantum-critical regime is relevant for 
electron-doped cuprates, in which antiferromagnetism 
emerges near the electron density at which the Fermi 
surface passes through (7r/2,7r/2), i.e , along zone di- 
agonal 2fci? — (tt, tt) (Refs.[22, l23(). We computed the 
velocity renormalization in the QC regime and found 
uj = i,,{l- 0.82(efc/T)i/4), where efc ex (fc - kp)*^^. Ob- 



serve that a given k, ut still increases with increasing T. 

To conclude, in this paper we considered temperature 
dependence of the nodal Fermi velocity in 2D systems, 
vf{T) = vf{0){1 + S{T)) We have found that S{T) is 
positive and linear in T in any Fermi liquid. The linear 
in T term comes from the screening of the interaction 
by just one particle-hole bubble. The corrections due 
to extra dynamical screening only give rise to regular, 
T^ terms. The slope of the linear term is quite large in 
the cuprates, and agrees with the measurements on opti- 
mally doped Bi2212. In the quantum-critical regime, the 
correction to the zero-temperature dispersion is singular, 
and scales as l/r°i'^ at the lowest T, and as l/T°-'^^ at 
intermediate temperatures. 

The experiments by Plumb et al [9| have been per- 
formed only at optimal doping. In the theory, the prefac- 
tor of the T term increases with decreasing doping, hence 
the prefactor should decrease in the overdoped regime, 
and increase in the underdoped regime. The measure- 
ments of the doping dependence of the slope are clearly 
called for. 
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